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Abstract 

Employing a recently developed technique of component weighted two parti- 
cle Green's functions in the CPA of a binary substitutional alloy A C B\- C we 
extend the existing theory of excitons in such media using a contact potential 
model for the interaction between electrons and holes to an approximation 
which interpolates correctly between the limits of weak and strong disorder. 
With our approach we are also able to treat the case where the contact interac- 
tion between carriers varies between sites of different types, thus introducing 
further disorder into the system. Based on this approach we study numeri- 
cally how the formation of exciton bound states changes as the strengths of 
the contact potentials associated with either of the two site types are varied 
through a large range of parameter values. 

PACS numbers: 78.50.-w 78.20.-e 71.55.-i 71.10.-x 
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I. INTRODUCTION 



Excitonic optical absorption in strongly disordered semiconductor alloys has been of great 
interest to both semiconductor physics and technological applications of semiconductors and 
semiconductor structures in the past. Although the theory of an optically excited electron- 
hole pair scattered or bound under the influence of the mutual Coulomb interaction into 
an exciton is quite well understood in weakly disordered systems and even at finite carrier 
densities in the respective conduction and valence bands, theories for strongly disordered 
systems such as alloys of insulators and some semiconductors are still very incomplete. 

One of the few theories of disorder which interpolates correctly between all regimes of 
disorder strengths and concentrations is the coherent potential approximation CPA first de- 
veloped by SovenS and Taylor! and subsequently extended to a proper two particle theory 
by Velickyi. The CPA predicts correctly the occurrence of split-off impurity bound states 
and whole impurity bands, once the relative disorder strength becomes of the order of the 
width of the unperturbed single particle density of states considered. Despite the numerous 
advantages that the CPA provides, it has proved difficult to incorporate a treatment of a 
carrier-carrier interaction into its framework and to our knowledge the only attempt to find 
a joint treatment of both effects has been made by Kanehisa and Elliott! who introduced 
a random phase like decoupling of the disorder average in the corresponding perturbation 
expansion in combination with a contact potential model for the Coulomb interaction. Al- 
though this description produces a number of correct features in a limit of low disorder, it 
could not be extrapolated successfully to the case where the disorder becomes stronger, i.e. 
when the joint density of states splits into two components and the excitons move to a more 
localized Frenkel limit, where a separate approach is needed. No theory has been available 
for the region of intermediate disorder strengths where the bands are about to split. 

An improved treatment of this problem has become possible through the development 
of a properly weighted two particle CPA in a previous paper! which allows to distinguish 
between different site types involved in the propagation of two particles during an absorption 
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process. We will show that differentiating between alloy components in a perturbation 
expansion with an appropriate decoupling procedure of the disorder average can yield the 
correct interpolation to all ranges of disorder for a contact potential model. Furthermore, 
this allows us to distinguish different contact interaction strengths between the carriers if 
these meet on different types of sites, whereby these strengths are modeled to represent the 
electron-hole interaction within an atomic radius. 



II. ABSORPTION IN THE NON-INTERACTING DISORDERED ALLOY 

A. CPA Model of the Disordered System 

We summarize here the results for a component weighted two particle CPA. Our model 
system is a binary substitutional alloy of components A and B with respective concentrations 
c and 1 — c, on a simple three dimensional monoatomic lattice. Both components are assumed 
to exhibit a direct band gap in their pure phase. Our two-band model Hamiltonian without 
carrier-carrier interactions can thus be written in a site representation as 

H = J2iW4+lCn + W?d + n+l d n } + Y.{U C nC + n Cn + U V n dtd n } (1) 
l,n n 

where the c n and d n are the annihilation operators for an electron and a hole on the site 
n, respectively. W 11 is the periodic part of the Hamiltonian transferring particles between 
different sites in the band fi and U£ is the matrix element in the respective bands which 
assumes the values U% G {£%£ > b\ depending on whether n is an A— or a I?-site. The 
disorder in the system is generated by the difference of the on-site energies V 1 = e^ A — e^ B in 
the respective bands. 

The CPA for a disordered medium is introduced by the usual methodBi'i of placing the 
impurities in a self-consistent medium such that the scattering off a single impurity vanishes 
on the average. Independent of this approximation, the single particle propagator of the 
disordered medium G relates to the one of a pure medium g through the Dyson equation 

G M = g" + gWG». (2) 
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The the self consistent CPA condition requires that the average propagator G M of the effective 
medium fulfill the relation 



G» = g» + g^G». (3) 

where E^ 1 is the CPA single particle self energy of the effective medium. S M itself is de- 
termined through the single site condition that the average atomic T-matrix of the self 
consistent medium Tj" defined as 



T% = (U£ - W) + (U£ - (4) 

be zero. Here we have introduced the site diagonal single particle function F = G(n,n). S M 
and therefore satisfy the self-consistent relation 

n = m\ ±1 ~ SM ] + (±-c)le B -m 

Weights can now be attributed to the single particle functions before averaging through 
applying operators S*/ h to them which effectively exclude either impurity or host sites from 
the choices of start or end sites of the propagation of the particles. We make the following 
definitions 

\ttii _M1 /T/tt. r<A rzAr<- r~<AA wAri "A. riAB 'aAri "B 

i n Bsl i ° —"a Lj - U — "itW^t' U 1— 'u M 1— 'ti ' fr% , 

(6) 

_ rriii i\ru. nB _ ^Br<- nBB _ ^Br< wB. c<ba _ ^Bn ^a 

— FA ~ U n\l V J ^ — "At^J ^ — ^n^fi^ni " — -^A*"/*) 

Averaged versions of such weighted single particle functions have been introduced by Aiyer et 
al.@ It will be shown that they are particularly useful in their unaveraged form when employed 
in the definitions of two particle functions in order to calculate weighted components of the 
polarizability which can thereupon be used in the treatment of the excitonic absorption. 



B. The Average Polarizability in Linear Response Theory 



In the following we consider a linear response expression for the polarizability of the Kubo 
type. The particular kind of two particle function needed in this case is determined by the 
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form of the dipole operators which account for the interaction of the electronic excitation 
with the radiation field and which in the case of allowed interband transitions! effectively 
couple the two single particle resolvents in the respective bands: 

= ^2 I nc)p(nv | . (7) 

n 

The total optical polarizability in a Matsubara representation which we adopt for conve- 
nience can be written as 

(rioM) = -/T 1 ^(Tr v {i& c {iz)&G v (iu - iz)}), (8) 

iz 

where the angular brackets denote the configurational average, (3 is the usual inverse tem- 
perature and Tr v denotes the trace over the band v as in Ref. |9| 

Tr v (- ■ ■) =Y J (nv \ ■ ■ ■ \ nv). (9) 

n 

Here it is assumed that the dipole matrix elements are taken for the transition from a valence 
band of p-wave symmetry to a conduction band of s-wave symmetry in which case they can 
be assumed to be essentially constant and we normalize them to unity This implies that for 
the calculation of a properly averaged polarizability which is void of any further interactions 
we need to consider the following two particle function 

K(zi,z 2 ) = ^2K(z 1 ,z 2 ;n,m) = ^2(nc \ (G c (zi) \ mc){mv \ G v (z 2 )) \ nv). (10) 

m m 

Apart from the cumulative function K stated above it is now also possible to consider its 
weighted components by using the appropriately weighted versions of the single particle 
functions in their definitions. 

K and its weighted forms are obtained by means of the two particle extension of the CPA 
by VelickyS and our method of obtaining the appropriate weig The concept is to use the 
T-matrix representation of the unaveraged unweighted or weighted single particle functions 
from (|6|) to obtain a CPA vertex correction A(zi, z 2 ) and average weights for the correlated 
motion of two particles in terms of averages over products of atomic T-matrices. With the 
vertex correction, the equation for the correlated two particle motion can be decoupled as 
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=1-A(*„ *)«(*,*)• (U) 
where -R(zi, 2^) is the average-decoupled two particle function 

R(zi, z 2 ) = G c (zi] n, m)G v (z 2 ] m, n) (12) 

m 

and the vertex correction is found to be 

A(zuZ2) - V c F c { Zl )-V v F v { Z2 y (13) 

If one assumes that the conduction and valence band dispersions are similar in shape 

£ c(k) _ _^_ e v (k) 
w c w v 

one finds a more explicit representation for R(zi, z 2 )Q 

R(z z ) - w c F c{z\) ± w v F v (z 2 ) 

U 2 w c [z 2 -T, v (z 2 )]±w v [z 1 -T, c (zi)}' 1 

In Ref. [| we found that for the class of two particle functions as represented by K(zi, z 2 ) two 
types of single weighted functions and three different double weighted ones exist. Similar 
to the single particle theory it was found that these weighted two particle functions can be 
expressed in terms of the unweighted one if energy dependent weighting factors added to 
and multiplied onto it. 
Introducing 

,, s (1 - c)A(zi, z 2 ) 

azu Z2) ~ (K-s^OXK-s^))' (16) 

cA{z u z 2 ) 

' ( *'* )= E.WE.W <17) 
= AM *' l)FM , (18) 



we find that 



K A = £K; K B = rjK, 

R AA = £2 R + T; R BB = V 2 R + 7; R AB = R BA = ^ R _ ^_ 



(19) 
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These functions can now readily be used to obtain corresponding components of the polar- 
izability which are calculated the same way as in @ 

(n* y (^)> = -/r 1 Y, k xy (iz, iuj - iz), (20) 

iz 

where X and Y denote possible weights (X, Y G {A, 5,0}). Functions of this type have 
been discussed analytically and numerically in Ref. [5] for various regimes of disorder. 

III. INTRODUCTION OF ELECTRON-HOLE INTERACTION 

A. Difficulties with the Treatment of the Coulomb Interaction in Disordered 

Systems 

The inclusion of a carrier-carrier interaction into a model for the optical polarizability is 
essential for the treatment of exciton effects. However, an analytic treatment of this problem 
in strongly disordered systems such as alloys of insulators and of semiconductors using the 
true or screened Coulomb interaction seems to be almost impossible. This is due to the fact 
that the corresponding Bethe-Salpeter equation 

(}C(iuj;i,j;k,l)) — (JCo(iu;i,j;k,l)) + (JCo(iuj;i,j;a, b)u(a, b;c, d)JC(iuj;c, d; k, I)). (21) 

a,b,c,d 

is extremely hard to decouple. Here u(a, b; c, d) = v(a—b)6 ac 5bd— , w(a—c)S a bS c d is the Coulomb 
interaction in the site representation with its direct and exchange part and JCo(iu;i, j; k,l) 
is the non-interacting configuration dependent two particle function 

JC (iLu;i,j;a,b) = -/T 1 ^ | G c {iz) | k)(l | G v {iz - iu) \ j). (22) 

The exactly averaged polarizability including exciton effects is calculated from the interact- 
ing version of this function (^Tj) as 

(n(iw)) = X)Ww;i,z;ife,A;)). (23) 

i,k 

The main difficulty in solving equation ([H]) arises from the fact that the disorder average 
over the second term creates higher order correlations in the particle motion than those 
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which are accounted for by the inclusion of the vertex corrections. These additional average 
induced correlations lead to an effective Coulomb interaction which is renormalized in a very 
complicated fashion. 

For these reasons two drastic simplifications have been made in earlier treatments of the 
problem. First the Coulomb interaction is replaced by a short range contact potential which 
binds at most one state in a joint band limit when the disorder is relatively weak. This state 
is meant to represent the 1-s exciton from the Coulomb series which usually dominates the 
optical absorption spectrum below the continuum absorption edge. Secondly the disorder 
average in (|2lD is decoupled into a product of averages. The average polarizability including 
this carrier-carrier interaction is then obtained as 



where u is the uniform strength of the contact potential and the average was replaced 
through a horizontal bar over corresponding quantities. However, this type of decoupling of 
the average becomes reasonable only in the limit of weak disorder, when the electron-hole 
pair experiences many impurities before its constituents are Coulomb scattered (Wannier 
excitons). Effectively, therefore no higher than two particle correlations can be accounted 
for within this approximation and it will be shown that for stronger disorder even these 
correlations are significantly misweighted and overcounted. It hence must clearly fail in a 
regime of strong disorder where excitons may be primarily associated with two separate 
impurity bands. 

For the strong disorder limit it is shown explicitly in Ref. [| how the solution for the two 
particle function is modified towards a regime of localized excitations (Frenkel excitons), 
which is equivalent to letting the overlap of the atomic wavefunctions go to zero, while at 
the same time the disorder strengths are kept constant. If fl2~l~|) is considered before averaging 
in this limit, it goes over to 

[1 + 7r (ico>; i)v (0)]/C(iu;; ii, kk) = dik^oiiuj] i) + vr (iu;; i) ^^w{i — j)JC(iu; jj, kk), (25) 




(24) 
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where no(iu;i) is the atomic polarizability on the site i and w(i — j) is the exchange part 
of the Coulomb interaction which is now entirely responsible for electron transfer from one 
atom to another. It is this term which makes an exact solution of equation ( p5j ) still difficult. 
If now a the exchange w is assumed to go to zero as well we obtain the limit of isolated 
atoms. 

The solution of ( 25) then becomes 

JC(icj; ii, kk) = (26) 

which can be averaged exactly and summed to yield 

cn (iiJ,A) (1 — c)7i (iu, B) 



(27) 



l + v(0)7r (iuj,A) l + v(0)ir o (iu,B) 

where A and B within ttq denote here that the atomic polarizabilities are considered on 
either an A or a B atom. Since the present limit corresponds to letting the bandwith of 
the components of the joint density of states go to zero, given some arbitrary strength of 
disorder, the same limiting behavior of the polarizability can be attained! by letting the 
disorder strengths become much larger than the band width which yields two energetically 
separated contributions to the joint DOS and the formal properties of fl2"T| ) are recovered. 

The appearance of (|27D suggests that one should also allow for the possibility of using two 
different contact interaction strengths on different types of sites thus allowing for different 
types of atomic screening in different alloy components. 

It is clear that equation (|27|) is not obtained by the extrapolation to a strong disorder limit 



from equation (|24| ) and so far there had been no indication as to what type of approximation 
one should use in the description of an intermediate disorder regime, where the strengths of 
the disorder are of the order of the corresponding half-widths of the single particle bands 
w M , which at the same time extrapolates correctly to the asymptotic behavior predicted by 
equations ( p4|) and ( |27| ) for weak and strong disorder, respectively. 
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B. Component Weighted Approximation of the Scattering Expansion 



Based on our results in Ref. || we employ in the following a component weighted scattering 
expansion to obtain the excitonic polarizability, while at the same time we allow the contact 
potential strengths to assume different values on different types of sites. We use a similar 
decoupling of the disorder averages to obtain a new approximation for the scattering in 
an intermediate disorder regime. Subsequently, we show that this approximation renders 
exactly the interpolation behavior that we had hoped to achieve. 

We introduce the following attractive weighted contact scattering interaction strengths 



ua 
ua = — ; 

c 



u B 



u b 



(28) 



associated with the A— and .^-components of the medium, whereby ua and ub represent 
bare scattering strengths, and the limit of uniform interaction strengths corresponds to ua — ► 
Ub — > u. If one distinguishes now the different types of sites involved in the propagation of 
the two particles during a scattering process, the expansion of the scattering series can be 
represented pictorially as in Fig. |I|. Mathematically this corresponds to writing 



+ 



A~ t\AB~ 



iB 



B~ TtBB- 



iB 



(29) 



where both, the diagrams in Fig. [I] and the terms written down in (^9]) above represent the 
expansion to second order. These and higher order terms can can be conveniently rewritten 
in a 2 x 2 matrix scattering formalism by introducing 



(30) 









( ~A 
-U A 


^ 




f f\ AA 


tjAB \ 


n = 


u = 






n = 




V n ° J 




K ° 


-u B J 




frAB 


fvBB 



Equation (|29|) then goes over to 

n = n + n^urio + n^un uno + 



(31) 
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The matrix products can be summed to give 

oo 

n = n + n[u^ [n u 



n=0 



n , 



(32) 



which can be evaluated as a matrix geometric series as 



CXJ 

E[n u 



?1=0 



1-IInU 



-1 



(33) 



and equation fl3~2|) is calculated to yield 



tt _ frA i ttB u A (Ujni + ggnn + u B (TlE?(l + u A U AA ) - 2U A u A U AB u B U B 

0+0 (i+u A nr)(i+u B n* B )-u A u B (^ B ) 2 ' 1 j 

Equation (0) is the central result of this work and much of the further discussion will be 
based on its properties. 



C. Limiting Behavior 

In order to test the usefulness of equation fl3"4|) it is necessary to investigate its behavior 
in various limits of disorder. 



1. Weak disorder limit; uniform interaction 

To establish a connection with formula (^) we first study the limit of low disorder and 
asymptotically equal bare contact interaction strengths u A = u B = u. It is well knownB 
that in this limit the CPA self energy goes over to the virtual crystal limit S c /„ — > cV c / v 
and the vertex correction to A — ► c(l — c)V c V v . Accordingly, the weights in the two particle 
function defined in (|16|)-(|T8|) go over to £ — >• c, r\ — >• (1 — c) and 7 — > c(l — c)F c F v . As a 
result of this it is possible to pull the two weights £ and rj multiplied with the unweighted 
two particle function out of the energy convolution, since they are now independent of 
energy. The impurity weighted versions of the polarizability can therefore be expressed 
in terms of the unweighted function as Tl A = cfl and Yi AA = c 2 fl + c(l — c)Q, where 
Q(iu) = — J^iz F c (iu)F v (iz — iu). The behavior of the host weighted functions and the 
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mixed one follow in analogy. Through inserting into (0) it is found that the scattering term 
factorizes as 

fi oU n (i + un) 

n n ° " (i + M n )(i + ^)' (35) 



which is identical to (|24]) . 



2. Strong disorder or split band limit 

In Ref. |5] we had shown explicitly that in the regime of strong disorder, when the widths 
of the single particle bands become negligible compared to the disorder strengths involved, 
the mixed component of the CPA two particle function K AB go to zero whereas the double 
weighted functions K AA I BB become effectively identical to the single weighted ones K A I B . 
The same behavior also translates to the weighted polarizabilities and therefore in this limit 
([34]) goes over to 

n A Ti B 

n = TT#w + r+^W' (36) 

Even though this already looks very similar to ([27]) it is not trivially the same. However we 
showed in Ref. [| that in this limit 

K >x ^^Cr c ^ crys (k)G v ^ crys (—k) — x^ K c J ys (37) 
k 

where x A l B is the concentration c and (1 — c) of the A— and 5-component, respectively, and 
the subscript crys is chosen to label the corresponding two particle functions of the pure A— 
and 5-media. This in connection with the definition of u A and u B from (]28| ) finally shows 
that ( p?P and (|34] ) are indeed identical in this limit. 

3. General behavior 



In order to examine the predictions that (131) makes for a general case, we introduce the 
following definitions 
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(38) 



where T(iu) = —(3 1 ^^^{iz^iuj — iz). Due to probability conservation the other weighted 
functions are defined as 



n* = (i- x )n ; = x (i - c)n - r ; n£ 



BB 



;i-2 X + x c)n + r. 



(39) 



In a general case, the polarization weights x an d C are only very indirectly related to the 
original weights £ and rj for the unintegrated two particle function K. They may assume dif- 
ferent values for any particular energy at which the polarizability is considered and therefore 
they generally depend on energy. Inserting this into fl34|) we find that 



+ u A u B 



n = n -ul{[ x 2 u A + (l - xfu 

x {l + [xCu A + (1 - 2 X + xQu B + u A u B T 



x(C-x)n + rj} 

n + {u A + u B )T + u A u B X (C - xMY 1 ■ (40) 



If this this is compared to a corresponding Dyson equation for the polarizability with a 
renormalized self-energy- like expression M 



IT = ITo + IToMn 



(41) 



this expression M, which can also be viewed as a disorder average "dressed" interaction, is 
found to be 

u A X 2 + (1 - x) 2 u B +u A u B x(l—x)® 



M 



(42) 



where we had introduced $ 



l + [u A + u B } X (l-x)$ 
x{( ~ x)So + r / [x(l — x)]- Introducing a generalized 



energy dependent average interaction strength 



u = X 2 u A + (1 - xfu 



2„~.B 



and a generalized contact potential fluctuation strength 



(43) 



A = X u A - (1 - X)u 



r.B 



(44) 
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equation (|42|) can be rewritten as 

M= -" + l + [ fi + A(l-2x)]» - <45) 
Finally introducing the renormalized function \1/ = $/(l + tt$) one obtains 

M — — u + X(1 A T X)A2tf T . (46) 
l + A(l-2x)« 

We argue that this result represents a generalized form of an average T-matrix approximation 
(ATA) of the polarizability of the medium with respect to the disordered contact potential 
relative to its "mean" background- value u which is taken as the basis of an energy dependent 
background medium similar to a virtual crystal (VCA) energy. To motivate our comparison 
more, we recapitulate the main features of the ATA for a single particle propagator G 
in a binary disordered medium with on-site potentials — e A and — e B on the A— and B- 
components, respectively The contact potentials have been defined with a negative sign 
here to match the definition of u A and u B from before, which were introduced as attractive 
interactions on an equivalent footing. The averaged T-matrix equation for G can be written 
as 

G = G + G f G , (47) 
in which G is the virtual crystal propagator 

S ° = TTT 9 - (48) 

where we have set e = ce A + (1 — c)e B and g is the the propagator for a medium without 
the on-site potentials. The average T-matrix T is defined as 

l + lc A - i]F a + 1 + [e" - e]F„ ' 1 J > 

where F = Go(m,m). Introducing furthermore 5 = e A — e B , the corresponding ATA 
selfenergy Y^ata relating g and G can be written as 

T c(l-c)5 2 Fo 
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The aforementioned analogy hence builds on the correspondence of the quantities 

e < — ► u 5 < — ► A 

g < ► Ilo ^ATA *~ -> M (51) 

G < > A F < > iff 

This correspondence becomes exact in the asymptotic limits discussed before and continues 
to hold qualitatively in an intermediate regime. The difficulty in finding a rigorous compari- 
son in the most general case stems again from difficulty of finding general expressions for the 
energy dependence of the integral weights x an d C an d the integrated diagonal correction 
term T. 

The reason why our weighted expansion scattering expansion can be expected to yield 
a better result than an undifferentiated decoupling of the Dyson equation for the contact 
potential, as it is shown to do by our numerical calculations in the next section, can be 
understood from the argument that averages over higher moments of n which occur in 
this expansion such as (1T5) and (IIq) are better approximated by a decoupling as (ITq) ~ 

(n£)7c+(n^) 2 /(i-c) an d (ng) ~ (n^) 2 n^/ c 2 +2n^n^/c(i-c)+n^(n^) 2 /(i- c ) 2 

rather than (IIq) ~ IIq and (IIq) ~ IIq, respectively, because the former approximation re- 
duces a wrong weighting and overcounting of scattering processes that the latter approxi- 
mation erroneously infers. 

IV. NUMERICAL RESULTS 

In this section we discuss a numerical implementation of our previous results within a 
commonly used model in three dimensional systems. Since within the quasi-ATA contact 
potential model of the electron-hole interaction the excitonic polarizability can be obtained 
directly by inserting the weighted components of the average polarizability for the free 
particle system into equation ([34]) , we draw strongly from the results obtained for the non- 
interacting system in Ref. [5[ 
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The main input necessary for the numerical implementation of a single site CPA are the 
single particle densities of states for the conduction and valence band. As in Ref. |5] we take 

(52) 

p^(E) = \E\< Wfl . 

It is well known that this density of states does not exhibit any of the finer structure of 
a real system, but it includes the more global features of a wide class of systems giving a 
finite band width and the appropriate van Hove singularities at the band edges. Moreover, it 
matches well with the spirit of a single site CPA which, even though it reproduces the global 
influence of disorder on the system correctly in all regimes, cannot account for more detailed 
structure brought about by the scattering of particles off clusters of impurities which usually 
involve momentum dependent self energy contributions whose evaluation would require a 
more specific knowledge of the single particle dispersion laws. 

For amalgamation type solid solutions, i.e. alloys whose densities of states of the pure 
A— and 5-substances overlap to a great extent in the respective conduction and valence 
bands, CPA results already exist^ for the case ua = ub = u. We shall therefore focus on the 
region of parameters which corresponds to intermediate disorder, where our approximation 
provides significantly better results than earlier approaches. A restriction to a rather specific 
region of parameter space becomes necessary since we now have six independent parameters 
which govern the behavior of our spectra, i.e. the concentration c, the relative valence - 
conduction band width w v /w c , the disorder strengths V c and V v in the respective bands and 
the corresponding disordered contact interaction strengths ua and Ub, accounting for the 
dynamic part of the electron-hole correlations. 

In determining the region of interest we note that, as pointed out by many workers 
(cf. for example the review article by Rashba0), the CPA produces a "pseudo" gap in the 
single particle DOS at all concentrations once the disorder strengths exceed the half-width 
of the single particle bands involved, > w^, and also translates to the joint DOS in some 
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regimes. In more realistic systems, however, usually a true gap only exists once the disorder 
strengths exceed the full band widths > 2w^, even though in the regime 2w^ > > 
the number of states in the region of the pseudo gap is strongly suppressed. Because of 
this, we shall concentrate our attention to the beginning region of the true gap behavior and 
choose | |= 3w M . 

In calculating the absorption spectra, we have first obtained the weighted components 
of the polarizability without electron- hole interactions. Figures |2| (a), (b) and |] (a), (b) 
show the real and negative imaginary part of Ho and all its decompositions into single 
and double weighted components for parallel disorder sgn(V c ) = sgn(V v ) and anti-parallel 
disorder sgn(V c ) = —sgn(V v ), respectively. Specifically, the conduction and valence band 
disorder have been chosen to be V c = 3.0 and V v = ±2.4 in connection with the half- 
bandwidths w c = 1.0 and w v = 0.8. It should be noted, that all energies occurring in our 
results effectively scale with the half-width of the conduction band w c as it is normalized to 
unity. 

We can see from Fig. |2] (b) and Fig. |] (b), that in the region of parameters considered, 
the joint density of states has already split quite clearly into distinct A— and 5-parts for 
both cases of disorder and the difference of the double and single weighted components has 
become relatively small due to the dominance of the total diagonal matrix elements of the 
two particle function in this regime. Note that in all the figures for the non-interacting 
polarizability only the single ^45-function has been plotted and the components of the 
spectra sum as Ilo /B = Uq A/bb + fi£ B and fl = fi£ A + Tl BB + 2Tl£ B . 

In the case of parallel disorder of Fig. Q (b), the joint DOS is divided up into into a central 
bulk part which is mainly constituted from AA— and 55-transitions and a separate set of 
flanks on either side of it, which contain a notable amount of AB— transitions. In the case of 
anti-parallel disorder of Fig. [| (b), on the other hand, it decomposes into two separate parts 
with a pronounced gap between them, whereby the upper part consists mainly of AA— and 
the lower part of 55-transitions. Even though the joint function appears to be very close 
to zero in the gap region for this case, we discover a finite contribution of double weighted 
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components centered between the split bands!. These states can be considered "mute" in 
the absence of any electron-hole interaction since cumulatively they do not contribute to the 
polarization. 

It appears from these results as if the system is already very close to the strong disorder 
asymptotic behavior where the split components are essentially independent, but we shall 
find that this is indeed not the case, agreeing with the result of Onodera and Toyozawa0 
who state that this regime is reached only beyond a relative disorder strength of about 
V^/w^ = 10. If in the interacting case states are pulled down from the upper band into this 
region by means of the electron-hole interaction, they will experience a significant broadening 
and hence change the excitonic absorption in a large region as as will be apparent from the 
plots in Fig. |. 

In order to show the effect of the various components of the densities of states on the 
formation of excitons in a concise way we have plotted stacks of three dimensional overlays 
of the obtained spectra. In every separate stack of the spectra in Figs. |] (a)-(d) and Figs. || 
(a)-(d) the strength ub of the contact interaction on a -B-site of the alloy - here associated 
with the lower lying component - is kept constant whereby the interaction strength ua on 
an A-site - associated with the higher lying component - is varied through an interval of 
strengths from zero to a value where it is sufficiently large to pull the exciton below the lowest 
component of the contributing unperturbed joint density of states. The same procedure is 
performed over all stacks of plots in Figs. || and |5] as ub is gradually increased through a 
similar interval of values as ua- 

Although, in an experimental situation it would probably be much easier to vary such 
parameters of the system as the concentration and to some extent also the strengths of 
the disorder by using different substances for the production of the solid solutions, a strong 
variation of the contact potential strengths ua and ub is better suited for a theoretical study 
of the global features of the excitonic absorption predicted by our model. As a result we 
are able to analyze the excitonic absorption through a whole region of the three dimensional 
parameter space spanned by u, ua and ub- In addition, since it is difficult to find a view- 



angle of the resulting plots which shows all features of the spectra simultaneously, we have 
cut off the resonance peaks at an appropriate value and overlayed a contour representation 
of the spectra which renders more detailed information about the position and the width of 
the excitonic resonances as the interaction strengths are increased. 

As can already be seen from equation Q53|), the width of the resonances will be largely 
determined through the behavior of the double weighted components, i.e. through the 
magnitude of the imaginary parts of TLq A and TIq B as well as fl^ B and at the solutions of 
Re[n^ A ] = —1/ua and of Re[n^ B ] = —1/ub, respectively. This implies that it is possible to 
have a very sharp resonance of the exciton peak associated with the higher lying component 
of the joint DOS deep within the region where the DOS is relatively large, but almost entirely 
consists of the opposite component. On the other hand the resonance can be very broad 
even in a region where the cumulative DOS is practically zero, due to the presence of finite 
double weighted components which mutually cancel out to a great extent, once they are 
summed. 

In Fig. |] (a) the set of spectra commences with the completely unperturbed joint DOS 
for parallel disorder in Fig. |2|. The A-exciton which is being pulled out at the lower end of 
the central bulk part, which consists almost entirely of A A— and 55-components, becomes 
very sharp as soon as it leaves the pure A-component at about u = —0.7 corresponding to an 
^-interaction strength ua = 0.45. Beyond this value almost the entire oscillator strength of 
the A-component is found in the resonance and thus only the pure 5-contribution remains 
in the central part as ua is increased further. It should be noted that these resonances show 
a finite width also when they are exterior of any of the DOS contributions due to a small 
artificial imaginary part of about 10~ 3 units which has been added to the energy in order to 
ensure the correct analyticity of the quantities involved at the given numerical accuracy. A 
major broadening occurs subsequently as the resonance crosses the lower AB-iegio\i which 
sharpens again once it has passed its lower end. The overlayed contour plot suggests that 
in the region without almost any AA-states, the increase of the excitonic binding energy is 
throughout linear with the increase to the interaction strength ua- However, as the bound 
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state passes through the the lower split off flank, which is constituted of about 50% AB— 
and 25% AA— and 55-components, respectively, - see 10 x enlarged region in Fig. (b) 
- it gets broadened and the propagation path of the bound state seems to attain a parallel 
shift (of about u = 0.7) with respect to the initial one, corresponding to a constant addition 
to the binding energy beyond the lowest flank of the spectrum. This additional deepening of 
binding is found through all stacks of spectra for parallel and anti-parallel disorder equally. 
It can be observed throughout Fig. || (b)-(d) that as the A-resonance passes the lower AB- 
flank most of the states from this region are being absorbed into the resonances while the 
same is true for the the upper AB-Q&nk of the spectrum which feels the effect from afar. A 
very similar behavior is observed as ub is increased over the series of stacked plots and the 
last bits of the upper flank get absorbed in this process. Passing through the lower AB-fiank 
also gives a strong broadening which accounts for the persistence of S-states which are not 
influenced by the A-interaction. 

The plots also show that the formation of excitonic states associated with either of the 
underlying components of the alloy is largely independent of correlations between the two 
interaction strengths ua and Ub, i.e. the formation of bound states associated with one of 
the two material components is seen to be as good as unaffected by a variation of interaction 
strength associated with the complementary component. The only slight correlation which 
comes into play is via the mixed Ai?-components to which both ua and ub couple weakly. 
It can be observed that as the resonance associated with the higher lying component passes 
through the shifted body of the lower lying one, both this resonance and the peak width of 
the absorption edge narrow in the case of parallel disorder Figs. ^| (b)-(d) and broaden in 
the case of anti-parallel disorder Figs. || (b)-(d). 

Figs. |^ (a)-(d) for anti-parallel disorder exhibit many similar properties to the ones 
for the parallel case, but they also display quite a few novel features. The A-resonance 
is strongly broadened as it is pulled into the center of the gap, which the non-interacting 
absorption for this case from Fig. |] (b) exhibits. This broadening stems from the AA— and 
SB-components which still prevail in this region and which are compensated by the AS-gain 
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contribution - see 5 x enlarged region of Fig. ^ (b) . It is very curious to observe however, that 
even after the resonance has passed to lower energies, there is a finite hump remaining in the 
gap center which is subsequently bleached as the 5-resonance also shifts downwards. The 
development of this hump can be followed through Figs. [| (a)-(d) if one looks through the 
trough that the A-resonance forms in the central region of the gap onto the residual spectrum 
visible at the back. It mainly consists of .B-B-states which remain in this region whereas 
both the AA— and A-B-components become largely withdrawn. Looking at the contour 
plots of this arrangement suggests that, in addition to the constant increase in binding of 
the 5-resonance as increases, which is comparable in magnitude to the one observed 
for parallel disorder as the central AB-region is passed, the constant of proportionality for 
the deepening seems to have increased corresponding to a faster linear deepening of the 
A-exciton binding with increasing ua- 

The spectra of Figs. |5| (a)-(d) are also particularly suitable to visualize that a weighted 
scattering expansion used in the context of an average decoupling and a contact potential 
model with two distinguishable interaction strengths renders a much more accurate predic- 
tion for the formation of the excitonic resonances than the unweighted one of equation (^3|) . 
Equation fl2"l|) produces a resonance whenever Re[n ] — —1/u and Im[n ] is small or zero. 
In a split band case as shown in Fig. f| the real part is seen to have a zero in the gap 
where Im[n ] — 0. This means that an increase of the interaction strength u would lead 
to a bound state being trapped between the split bands at this zero of Re[n ] in the limit 
11 — > oo. This behavior of course is clearly wrong, since one can expect that the binding 
energy of any occurring excitonic resonance deepens - in this case the binding energy of the 
resonance associated with the higher lying A-band - as the interaction strength increases. 
This failure in giving the right asymptotic description can be ascribed to a misweighting 
and overcounting even of scattering processes which contain two particle correlations only. 
Our approximation overcomes this problem to a great extent and it is void of the erroneous 
trapping of the resonance, which as expected passes on to lower energies and finally appears 
below the B-band as ua is increased, independently of the value of ub- 
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The variation range of the carrier interaction strengths - particularly the one for ua - 
which we have plotted throughout Figs. |3| and is very large and is unlikely to be found in 
semiconductor alloys. However, mixtures of molecular crystals which form extremely narrow 
bands and at the same time have rather strong band offsets might be candidates for such 
behavior. 



V. COMPARISON TO EARLIER WORK 

Much earlier work on excitons in strongly disordered binary solid solutions at finite 
concentrations focussed on systems where the exciton is tightly bound with a binding energy 
greater that the narrow band width. In this case the Frenkel exciton has been frequently 
described by a "single" particle Green's function based on the theory of molecular excitons 
by Davydovli!. One of the standard assumptions within this framework, which is based on 
the Davydov theory, is that the matrix elements in the optical absorption, m^k' = ( c , k | 
i ■ p | v, k'} are strongly localized in momentum space, i.e. = const x 5k,k'^k,o which 

means that the contributing transitions are not only required to be vertical due to the 
approximate absence of total momentum of the exciting photon but also to only occur at 
one single point in the bands at k = k' = (Davydov component), due to the momentum 
of the relative electron-hole motion being approximated to zero. 

Early treatments of the absorption of mixed molecular crystals employed the average 
amplitude approximation (AAA) first introduced by Broude and RashbaEl, which is able to 
roughly predict the position and gravity center strengths of the dominant transitions only. 
With the introduction of the CPA which has some important advantages over the AAA, 
such as being able to produce an approximation of the actual shape of the joint DOS, it 
became possible to implement the Davydov theory within a self consistent framework. 

This has been pursued by Onodera and ToyozawaS and a little later by Hong 
and Robinson0 who tried to model their results particularly on observations made in 
experimentsEUll on naphtalene and anthracene. The approach of all these authors, de- 
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spite its self consistency, results in serious restrictions to the predictions their model can 
render about aspects of genuine two particle behavior and it contains no explicit mention of 
an electron-hole interaction. It seems to give reasonable predictions on some of the features 
of Frenkel excitons in the strongly insulating substances as investigated in referencesS'S, 
where the calculated positions of the Davydov components are compared with the ones of 
the experimentally observed peaks in the absorption spectra. 

A significant disadvantage of this single particle exciton theory is that only one "joint 
density of states" bandwidth and disorder strength can be considered, rather than two inde- 
pendent values for both of these parameters, associated with the underlying conduction and 
valence bands. A behavior of the joint density of states as displayed for parallel conduction- 
valence band disorder which can lead to the filling of the central part of the joint DOS by the 
states arising from the AA— and 55-transitions even in a split band limit, as shown in Fig. 
||, could for example never be obtained in such a simple single particle picture. Additionally, 
the binding energy of the excitons does not relate to any electron-hole interaction, which in 
a real system determines their relative position with respect to the continuum absorption 
edge. 

The first work which is based on the footing of a proper two particle theory known 
to us is by Abe and Toyozawa@ who employed Velicky's two particle CPA! to calculate 
the absorption in a non-interacting system with Gaussian disorder. Kanehisa and Elliott! 
thereafter seem to be the first authors to consider both a two particle CPA for the disorder 
and a contact interaction as a model for an electron-hole interaction, but they were only able 
to obtain an asymptotic solution in the regime of weak disorder, using the average decoupling 
of the unweighted polarizability in (|4]). As they pointed out, their results can be regarded 
as applying to the amalgamation regime and a comparison was made with experiments on 
III-V alloys. 

One should note that the approaches of Refs. and the present work make the 

assumption that the optical matrix elements are constant over the single particle bands 
involved, even though only vertical transitions at k = k' are allowed, which is contrasting 
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the Davydov theory used in the earlier theoretical models of Refs. 0,0, where only the zero 
momentum component k = k' = is taken. In a more realistic situation, the appropriate 
matrix elements will exhibit a more general k-dependence, but it does not present any 
conceptual difficulty to modify our present results to model such CclSGS clS well, once the 
actual dispersion laws for the conduction and valence bands are known. 



VI. CONCLUSION 

The present work is the most extensive to date which treats the correlated motion of 
pairs of interacting particles in a disordered medium. It draws from the extension of the 
CPA to two particle propagation developed by VelickyH and its extension to weighted two 
particle Green's functions given in our earlier paper!. 

The model introduced here achieves a better but still approximate treatment of the 
interference between the effects of impurity scattering and direct two particle (electron- 
hole) interaction and it provides an interpolation scheme through all strengths of disorder 
as well as the possibility to include site dependent variations in the direct electron-hole 
interaction within an ATA-like framework. It is therefore most appropriate for systems 
where the disorder causes a small but distinct band splitting, where previous asymptotic 
models fail to work. 

There are many disordered systems where such interacting pairs of particles play an 
important physical role and we believe the method can be extended to such situations. 
Modeling the excitonic absorption provides a particularly obvious and parmount application. 
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FIGURES 

FIG. 1. Diagrammatic representation of processes to second order included in the weighted 
expansion of the excitonic polar izability. The electron hole pair can be excited on an A— or a 
-B-site and subsequently be scattered on either of such sites with the effective contact potential 
strengths ua and ub, respectively. Every connecting line between two sites denotes a factor of a 
non-interacting weighted polarizablity IIq-^ (X, Y € {A, B}), where X and Y are chosen to match 
the type of the start and the end site of an arrow and all conceivable paths occurring between the 
excitation and the recombination of the electron-hole pair are being summed over. 

FIG. 2. (a) Real and (b) negative imaginary parts of the non-interacting polarizability for 
parallel disorder at c=0.35. All energies scale with the conduction half-band width w c which has 
been set to unity, w c = 1, while w v is taken to be w v = 0.8. Since this spectrum is taken at disorder 
strengths V c = 3.0 and V v = 2.4, i.e. at the beginning of the true gap region of the CPA, the single 
and double weighted functions coincide to a great extent, and two flanks of the spectra which to 
a substantial degree contain mixed ^4-B-components have split off sideways (cf. the lower end of 
the joint DOS in (b) where the values were enlarged by a factor of 10). However, it becomes clear 
from the spectra in Fig. 3 that the residual interference of A— and S-absorption amounts to a 
significant influence on the excitonic absorption as the corresponding electron-hole interactions are 
introduced. 

FIG. 3. Excitonic absorption spectra created by means of equation (34) from the weighted 
components of the non-interacting spectrum in Fig. 2. In each of the plots (a), (b), (c), (d) 
the contact potential strength ub associated with 5-atoms is kept at a constant value while ua is 
increased from to 6 from the front to the back of the set of overlayed plots. The overlayed contour 
plots show the paths and broadenings of the excitons as the parameter strengths are varied. The 
contour lines are taken at level heights of 0.1, 0.5, 1.0, 1.5, respectively. Comparing the plots (a) 
ub = 0.0, (b) u b = 0.4, (c) ub = 1.5, (d) ub = 3.5. one can see that mutual interference effects 
of both interactions are relatively small. The main such interference effect is that ^-resonance 
attains some broadening once the yl-resonance has passed it, cf. Fig. 3 (c). 
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FIG. 4. (a) Real and (b) negative imaginary parts of the non-interacting polarizability for 
anti-parallel disorder and parameters c = 0.35, V c = 3.0, V v = —2.4, w c = 1.0 and w v = 0.8. The 
anti-parallel direction of the band offsets causes the imaginary part (joint DOS) in (b) to split into 
two distinct components mainly constituted by A— and ^-transitions, respectively, separated by 
a wide gap. Even though the cumulative function is very close to zero in the central region, the 
double weighted components are finite there (cf. the central region in (b) where the values were 
enlarged by a factor of 5) and hence they lead to a substantial broadening if an excitonic resonance 
is pulled out of the higher lying A-component into this region, as can be seen from the excitonic 
spectra in Fig. 5. 

FIG. 5. Excitonic absorption spectra based on the non-interacting spectrum in Fig. 4. While 
in each of the plots (a), (b), (c), (d) ua is varied from to 8 to pull the A-resonance below the 
onset of the lowest B states, ub is varied as (a) uj, = 0.0, (b) ub = 0.4, (c) ub = 0.8, (d) ub = 1.8. 
Opposite to the parallel case shown in Fig. 3 the interference of the two interactions amounts to 
a narrowing of the -B-type resonance as the ^4-resonance passes it, cf. Fig. 5 (c). Both here and 
for the cases shown in Fig. 3 this can be understood to originate from the mixed ^S-contribution 
into which both ua and ub couple weakly. Similar to Fig. 3 the contour lines are taken at level 
heights of 0.1, 0.4, 0.8, 1.2, respectively. Note also that an unweighted contact potential model 
with an uniform interaction u of equation (24) would in this case erroneously predict a trapping of 
an excitonic resonance in the gap region where -Re[fto] = at about uj = 0.8, as u — > oo. 
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